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ON THE DISTRIBUTION OF THE MAXIMUM VALUE OF THE 
CHARACTERISTIC POLYNOMIAL OF CUE RANDOM MATRICES 

Y. V. FYODOROV AND N. J. SIMM 


Abstract. Motivated by recently discovered relations between logarithmically correlated 
Gaussian processes and characteristic polynomials of large random N x N matrices H from 
the Gaussian Unitary Ensemble (GUE), we consider the problem of characterising the distri¬ 
bution of the global maximum of Dn{x) := — log | det(a:/ — iJ)| as TV —>• oo and x G (—1,1). 
We arrive at an explicit expression for the asymptotic probability density of the (appro¬ 
priately shifted) maximum by combining the rigorous Fisher-Hartwig asymptotics due to 
Krasovsky |35j with the heuristic freezing transition scenario for logarithmically correlated 
processes. Although the general idea behind the method is the same as for the earlier con¬ 
sidered case of the Circular Unitary Ensemble, the present GUE case poses new challenges. 
In particular we show how the conjectured self-duality in the freezing scenario plays the 
crucial role in our selection of the form of the maximum distribution. Einally, we demon¬ 
strate a good agreement of the found probability density with the results of direct numerical 
simulations of the maxima of Dn{x). 


1. Introduction. 

The space of all A x A Hermitian matrices H with probability density function 


( 1 . 1 ) 


P{H) oc exp(-2ATr(A2)) 


is known as the Gaussian Unitary Ensemble (or GUE) [H |3Hl 03] . Here and henceforth the 
variance is chosen to ensure that asymptotically for A —)■ cxo, the limiting mean density of 
the GUE eigenvalues is given by the Wigner semicircle law p{x) = (2/^)^! — supported 
in the interval x G [—1,1]. The characteristic polynomial Pn{x) = dei{xl — H) of the matrix 
H constitutes one of the most basic quantities of interest, encoding all eigenvalues of H 
through the roots of pn{x). As one varies the argument x over an interval containing many 
eigenvalues for a given realization of the ensemble, the value of the polynomial Pn{x) shows 
huge variations by the orders of magnitude for large A, see Figure [LT] for A = 50 and Figure 
Ofor A = 3000. 
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Figure 1.1. A plot of a single realization of \pj^{x)\e for N = 50. 

The global maximum is marked with a red circle. 
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Figure 1.2. A plot of a single realization of 2 log with 

N = 3000. The maximum value is marked with a red circle. 
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The purpose of this article is to describe the statistical properties of the highest peak 
displayed by the modulus of the GUE polynomial |pAr(x)|, namely the probability density 
for the maximum value attained by \p]s[{x)\ over the interval [—1,1] on the real line as 
N ^ oo. Our main result is the following 


Conjecture 1.1. Consider the random variable 


( 1 . 2 ) 


M := 


max 

a:e[—1,1] 


2 log \pn(,x) \ -2E(log \pn(,x)\) 


Then in the limit N ^ oo we have 

(1.3) = 2 log(iV) - ^ log(log(iV)) - (1 + oil))y + o(l) 


where y is a eontinuous random variable characterized by the two-sided Laplace transform of 
its probability density: 


(1.4) 


£( 62 ^") 


^ + 3)G(s + 7/2)^ 

’ G{s + 6)G{s + l) 


where r( 2 ;) and G(z) stand for the Euler gamma-function and the Barnes digamma-function, 
correspondingly. The normalization G can be evaluated explicitly as 

( 1 - 5 ) ^ = 

where A is the Glaisher-Kinkelin constant A = = 1.2824271291.... 


Remark 1.2. The product form of the Laplace transform fll.4p offers an interesting inter¬ 
pretation of the above results. Noting that r(l -|- s) is the moment generating function of a 
standard Gumbel random variable G, we can write 

( 1 . 6 ) y = G-^y’ 


where y' is an independent random variable with two-sided Laplace transform 


In the end of the paper we provide convincing numerical evidence that this Laplace transform 
does indeed dehne a unique random variable y' . This immediately implies that the probability 
density of y is the convolution of a Gumbel random variable with y' . Such a convolution 
structure is expected to appear universally when studying the extreme value statistics of 
logarithmically correlated Gaussian fields, see the discussion around and after Eq. fll.lip . 


In recent years, much interest has accumulated regarding the statistical behaviour of 
characteristic polynomials of various random matrices as a function of the spectral variable 
X. To a large extent this interest was stimulated by the established paradigm that many 
statistical properties of the Riemann zeta function along the critical line, that is C(l/2 + it), 
can be understood by comparison with analogous properties of the characteristic polynomials 
of random matrices [321 lai [121 m [291 ig. 

For invariant ensembles [3H1 SS] of self-adjoint matrices with real eigenvalues, statistical 
characteristics of Pn{,x) depend very essentially on the choice of scale spanned by the real 
variable x. From that end it is conventional to say that x spans the local (or microscopic) 
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scale if one considers intervals containing in the limit N ^ oo typically only a finite number 
of eigenvalues (the corresponding scale for GUE in fll.ip is of the order of 1/N). At such 
scales, standard objects of interest are correlation functions containing products and ratios of 
characteristic polynomials, which show determinantal/Pfaffian structures P 1251H51 151 IHl IB3] 
for Hermitian/real symmetric matrices and tend to universal limits at the local scale. Similar 
structures arise for properly dehned characteristic polynomials Pn{G) = det {l — U of 
circular ensembles (like CUE, COE, and CSE) [3H] of unitary random matrices U uniformly 
distributed with respect to the Haar measure on U{N) (and other classical groups) [121 El US] , 
whose properties on the local scale are indistinguishable from their Hermitian counterparts. 

Next, when x spans an interval containing in the limit N ^ oo typically of order of 
N eigenvalues one speaks of the global (or macroscopic) scale behaviour. At such a scale 
properties of Pn{x) display both universal and non-universal features, the latter depending on 
the ensemble chosen. The study of characteristic polynomials at such a scale was initiated in 
[30] where it was shown that the function Vn{0) = —2 log | det(l — U e“*^)|, with U belonging 
to the CUE, converges (in an appropriate sense) to a random Gaussian Fourier series of the 
form 

OO 

(1.8) V(e) = V ^ , 

where the coefficients Vn are independent standard complex Gaussian random variables, 
i.e. E{n„} = 0, = 0 and = 1. The covariance structure associated with 

such a process is given by W.{V{9i)V{ 62 )} = —21og|e*^i — long as 9i ^ 92 - Such a 

(generalized) random function V{9) is a representative of random processes known in the 
literature under the name of 1 /f noises, see [271 ES] for background discussion and further 
references. 

Recently the study of the global scale behaviour was extended to the GUE polynomial 
Pn{x) in [211 by using earlier insights from [31] and [35] . That work revealed again a structure 
analogous to that of fll.81) . though different in detail. Namely, it was shown that the natural 
limit of D]si{x) := — log \pp^{x) \ + E{log |pAr(a;)|} is given by the random Chebyshev-Fourier 
series 

OO ^ 

(1.9) F{x) = '^—^anTn{x), a:e(-l,l), 

1 vn 

with Tn{x) = cos(n arccos(a;)) being Chebyshev polynomials and real a„ being independent 
standard Gaussians. A quick computation shows that the covariance structure associated 
with the generalized process F{x) is given by an integral operator with kernel 

OO 

(1.10) E{F{x)F{y)} = Y, -Tn{x)T4y) = -- log(2|a; -y\), 

n=l 

as long as X ^ y. Such a limiting process F{x) is an example of an aperiodic l//-noise. 

Finally, one can consider an intermediate, or mesoscopic spectral scales, with intervals 
typically containing in the limit N ^ 00 the number of eigenvalues growing with N, but 
representing still a vanishingly small fraction of the total number N of all eigenvalues. The 
properties of the characteristic polynomials at such scales were again addressed in [21] where 
it was shown that for the GUE, that object gives rise to a particular (singular) instance of 
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the so-called fractional Brownian motion (fBm) [33 US] with the Hurst index hf = 0, again 
characterized by correlations logarithmic in the spectral parameter. 

The discussion above serves, in particular, the purpose of pointing to an intimate con¬ 
nection between Gaussian random processes with logarithmic correlations and the modulus 
of characteristic polynomials at global and mesoscopic scales. The relation is important 
as logarithmically correlated Gaussian (LGG) random processes and helds attract growing 
attention in Mathematical Physics and Probability and play an important role in problems 
of Quantum Gravity, Turbulence, and Financial Mathematics, see e.g. [18]. In particular, 
the periodic 1// noise fll.Sp emerged in constructions of conformally invariant planar ran¬ 
dom curves |1]. Among other things, the statistics of the global maximum of LGG helds 
attracted considerable attention, see [15] and references therein. Particularly relevant in the 
present context are the results of Ding, Roy and Zeitouni [15] on the maxima of regular¬ 
ized lattice versions of LGG helds which we discuss informally below. Let Vm = be the 
d—dimensional box of side length N with the left bottom corner located at the origin. A 
suitably normalized version of the logarithmically correlated Gaussian held is a collection 
of Gaussian variables : v E Vn with variance = 21ogA^ -|- f{v) and covariance 

structure 

N 

(1.11) E{0Ar^„,07v,n} = 21og+ r- -+g{u,v), foT U ^ V E Vn 

\u — v\ 

where ln+(t(;) = max(lnt(;,0) and both f{v) and g{u,v) are continuous bounded functions 
far enough from the boundary of Vat. Now set Mn = max^gy^^ (j)N,v and thn = VdlogN — 
^ log log N. The limiting law of is then expected, after an appropriate shift and 

rescaling, to be given by the Gumbel distribution with random shift: 

(1.12) P{y) = lim Prob{MN > ^at - y) = E (| , 

N^oo I J 

where the distribution of the random shift variable z depends on details of the behaviour of 
covariance fll.lip for \u — v\ ~ A^ and \u — v\ ~ 1, see the detailed discussion in [TS]. The 
random variable Z = is related to the so called derivative martingale associated with 

the LGG helds [I5] whose distribution is however not known. Recently it has been shown 
that the recentering term tutv in (I1.12p also holds for a randomized model of the Riemann 
zeta function [2], proved by revealing a special branching structure within the associated 
logarithmic correlations. 

We see that our coniecture 11.11 for the maximum of characteristic polynomial of large GUE 
matrices fully agrees with the predicted structure of the maximum of LGG in dimension 
d = 1. Note that the expression (I1.12p implies that the double-sided Laplace transform of 
the density p{y) = —-^P{y) for the (shifted) maximum y is related to the density p{z) of 
the random variable z as 

(1.13) E (e^Q = j p{y)e^y dy = T{s + 1) j p{z)e^^ dz = T{s + 1)E (e^Q 

which is in turn equivalent to the Gumbel convolution in eq. fll.6p . In fact our formula fll.Zp 
provides the explicit form of the distribution for the derivative martingale of our model, thus 
going considerably beyond the considerations of na. 
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From a quite different perspective, processes similar to fll.Sp and fll.9p appeared in the 
context of statistical mechanics of disordered systems when studying extreme values of ran¬ 
dom multifractal landscapes supporting spinglass-like thermodynamics [HI |25l [271 E]. The 
latter link is especially important in the context of the present paper. The idea that it is 
benehcial to look at |pAr(6*)| as a disordered landscape consisting of many peaks and dips, 
and to think of an associated statistical mechanics problem was put forward in [HI IH] . It 
allowed to get quite non-trivial analytical insights into statistics of the maximal value of the 
CUE polynomial sampled over the full circle 6 G [0, 27r], or over its mesoscopic sub-intervals. 
This was further used to conjecture the associated properties of the modulus of the Riemann 
zeta-function along the critical line, see some recent advances inpired by that line of research 
in [2]. Some relations between between CUE characteristic polynomials and logarithmically 
correlated processes (in the form of the so-called ’’multiplicative chaos” measures introduced 
by Kahane, see [H] for a review) was recently rigorously verihed in [IS]. The case of CUE 
polynomials however remained outstanding. 

It is our objective in this paper to provide two separate means of supporting Conjecture 
11.11 First, we will provide careful and explicit, albeit in part heuristic, analytical arguments. 
Although our technique is inspired by the approach of [23] it contains new nontrivial fea¬ 
tures necessary to overcome challenges arising from the non-uniform eigenvalue density p{x), 
reflecting absence of translational invariance for the CUE at the global spectral scale (note 
e.g. the non-trivial recentering in fll.2p b All this makes actual calculation for the CUE much 
more involved in comparison to the CUE and the limiting random variable u above appears 
to be more complicated than its CUE counterpart. Secondly, we will test our Conjecture 
with numerical experiments for matrices of size N = 3000 and around 250, 000 realizations. 
This is especially important as part of our analysis is based on very plausible but as yet not 
fully rigorous considerations. Finally, is natural to expect that the same distribution should 
be shared by the maximum modulus of characteristic polynomials for Hermitian random 
matrices with independent entries taken from the so-called Wigner ensembles, see [19] . 

Before giving the detail of our procedure in the next section we need to quote the following 
fundamental asymptotic result obtained by Krasovsky [35] which will be central for our 
considerations: 



(1.14) 


(1.15) 


where 



(1.16) 


and G{z) is the Barnes G-function. Differentiating with respect to a, we deduce that 


(1.17) 


E(21og |pj,(i)|) = lV(2i2 - 1 - 21og(2)) + 0(log(iV)/lV). 


where we used that C"(0) = 0. 
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The most salient feature of the asymptotics fll.lSp is the product of differences on the 
second line, which when rewritten in the form 


( 1 . 18 ) 


exp 


log 12- 0:^)1 , 


can be looked at as evidence of the limiting Gaussian process fll.9p with logarithmic covari¬ 
ance fll.lOp in the background. We will however stress that naively replacing the (shifted) 
log \pN{xj) \ with the corresponding 1// noise fll.Qp is not a valid approximation as the factors 
(1 — in fll.lSp do play an essential role in determining the extreme value statistics of 

\pN{,Xj)\. Let us finally note that had we suppressed the factors C{aj) the faithful descrip¬ 
tion of log\pN{xj)\ would be that of the regularized LCG process with covariance fll.lOp . 
the position-dependent variance 21nA^ -|- 2 In |\/1 — and the position-dependent mean 
N{2x‘^ - 1 -21og(2)). 

Acknowledgements: We are grateful for helpful comments from Ghristian Webb during 
the preparation of this manuscript. We acknowledge support from EPSRG grant EP/ J002763/1 

“Insights into Disordered Landscapes via Random Matrix Theory and Statistical Mechanics”. 


2. Statistical mechanics approach to the distribution of GUE 

CHARACTERISTIC POLYNOMIALS. 

Following the ideas of [23] we recast the problem of computing the value of the global 
maximum of |pAr(x)| (with an appropriate shift by the mean value) as a statistical mechanics 
problem characterized by the partition function 

N 

(2.1) Z^(/3) = y y ^ dx, /3 > 0, g > 0 

with the ’’potential” (j)Nix) = —2 (log |p 7 v(t)| — Elog |pAr(a:)|)), inverse temperature /3 > 0 
and /9-independent non-negative parameter q. Specifically, if we define the associated ’’free 
energy” as = —/3~^ log Zn{/3), then 

(2.2) lim J^(/9) = min 07v(x) = 2 max [log |pAr(x)| — Elog |pAr(x)|]. 

/3^0O 

Note that if compared to a similar partition function for the GUE case the main new 
feature in fl2.ip is the factor p{x)'^. Although naively the presence of such a factor may seem 
irrelevant when taking the limit /3 ^ oo, we will actually see that it plays a very important 
role in supporting our procedure of extracting the free energy for /3 exceeding some critical 
value. 

Now we aim to compute the integer moments of the partition function defined in (12.11) : 

(2.3) E(^‘(/3)) = 

In the limit N ^ oo the leading asymptotics of the above integral can be extracted by 
replacing the factor E \PN{xj)\'^^'^ with its asymptotics from fll.151) . In this way one 
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obtains 

(2.4) 


'iV 


l+/3^ 


-1 


) / n {2\xj-Xi\) ^^"dxi. 

J ,=l i<,<i<A: 


mYw)) - 

After changing variables Xj = 2yj — 1 the integral above assnmes the form 


.. dxi 


k 


(2.5) ^k(p^+,+i)-2y^Hk-i) jj , Y[ \y,-y^\-^f^"dy,...dyk 

_ 2^(/^^+Q+l)~2j5^fc(fc— 1 ) o ( Y^ + q Y^ + q _g 

fe 2 ’ 2 ’ ^ 

with the qnantity Sk{a,b, — 7 ) being the well-known Selberg integral 


( 2 . 6 ) 


(2.7) 


St(a, 6,-?) := 


X j X d 


dxi... dx}^ 


n lib 

-^[0,1]'= jJi 

r(a +1 - {] - i)7)r(& +1 - (j - i)7)r(i - jt) 


n 

l<2<_7</c 


k 

n 

i=i 


( 2 . 8 ) 


r(a + b + 2-{k + j- 2)7)r(l - 7 ) 
-Sk{a,b, - 7 ) 


r^(i -7) 

It is easy to see that the fonnd expression for the partition fnnction moments 

is well-defined for any 0 < 7 = < 1 and for an integer k satisfying 1 < fc < 7 


m 

-1 


To nnderstand how to deal with the case k > we recall that Krasovsky’s asymptotic 
formula fl 1.151) is valid only when all of the differences \xi — Xj \ remain finite when N ^ 00 , 
and should be replaced by a different expression when \xi — xf ~ N~^. One can check 
that the divergence of the integral for A; > l//9^ is due precisely to the fact that these near 
degeneracies become important. Relying on our experience with the corresponding situation 
for the CUE [23] case suggests that taking into account the correct short-scale cutoff cures 
the formal divergence, but changes the asymptotics of the moments E(^^) with N\ namely, 
these become of the order of for k > /3~‘^ whereas they are of the order of for 

k < Y~‘^- Such a change of behaviour will lead to a log-normal (far) tail in the distribution. 
Note that for the CUE case, the above behaviour conjectured in [23| was validated by recent 
rigorous calculation m- There is no doubt that the same mechanism is operational in 
present case as well and will be validated by extending the theory of m from Toeplitz to 
Hankel case. Actually, as argued in [ 2 T| the moments with k > play only a secondary 
role when addressing the question of extreme value statistics which is controlled exclusively 
by moments with 1 < fc < Our next goal is to use the latter integer moments for 

restoring the associated part of the probability density V{Zp) for the partition function. 
This will be achieved if we manage to find the distribution for a random variable zy whose 
positive integer moments are given by 


E(zJ) = Sk{a,b, - 7 ), 


a = b = 


q + Y^ 


i = Y^ 


(2.9) 
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Such a task actually requires finding a way to continue analytically those moments to com¬ 
plex k. Below we will arrive at the required continuation by exploiting a relatively simple 
heuristic procedure suggested in [25]. Note that in a series of insightful papers [3^H0lHnH2] 
Ostrovsky developed a rigorous mathematical procedure of the required continuation which 
provides an aposteriori justification of the results obtained via the heuristic approach. 


2.1. Analytical continuation of Selberg’s Integral. One starts with finding a recursion 
satisfied by S'fc(a,6,7) for integer k which is suitable for the continuation. By writing 
( 2 . 10 ) 


k-l 

Y[T{a + b + 2 -{k + j - 8)7) 
i=i 


r(2 + a + 6— (A; — 2)7) 11^=1 ^(2 + a + 6 — (fc+j — 2)7) 
r(2 + a + 6 - (2fc - 2)7) r(2 + a + 6 - (2A; - 3)7) 


one sees immediately that 
(2.11) 

E(z^) _ r(a + 1 - (A: - l)7)r(6 + 1 - (A; - l)7)r(l - k-f)T {2 + a + b - {k - 2)7) 

~ T {2 + a + b-{ 2 k- 2)7)r(2 + a + 6 - ( 2 A: - 8)7) 

It is convenient to introduce the moments Mp{s) of the random variable defined for any 
complex s as Mp{s) = . We then have E{z^p) = M (1 — k), E(z^“^) = M (2 — k) and 

after identifying s = 1 — A: the recursion fl2.1ip takes the form 

(2 121 ^ r(l + g + 75)r(l + b + 7g)r(l + (s - l)7)r(2 + g + & + (s + 1)7) 

+ 1) r(2 + g + 6 + 2s7)r(2 + g + 6 + (2s + 1)7) 

which is now assumed to be valid for any complex s. It is convenient to further use the 
duplication formula for the Gamma function : 

022: —1 

( 2 . 13 ) r{ 2 z) = ^^r(^)r(z +1/2) 

Vvr 


to get rid of the argument 2s in the denominator. Indeed, we have 

r(2 + g + 6 + 27s) = 21 +“+^+t' 2 T (1 + 7 s + (g + 6)/2)r(7s + {a + b + 3 )/ 2 )/V^ 
r(2 + g + 6 + 7 ( 2 s + 1)) 

= 2i+“+'>+N2^+i)r(l + 7 (s + 1/2) + (g + &)/ 2 )r( 7 (s + 1/2) + (g + 6 + 3)/2)/^ 
so that fl2.12p assumes the form 

(2 141 ^/ 3 (^) _ r(l + g + 7g)r(l + b + 7 g)r(l + 7(s - l))r(2 + g + 6 + (s + 1)7) 

+ 1) r(l + 7 (s + 1/2) + (g + 6)/2)r(7(s + 1/2) + (g + 6 + 3)/2) 

TT 1 

^ 22(i+a+fe)+(4s+i)7 r(l + (a + b )/2 + 7s)r((3 + g + b )/2 + 7s) 

Recalling that according to fl2.9p in our particular case g = 6= | + ^we now use the 

parameterisation g = gi + g2/9^, & = and /9-independent constants gi,g2,&i.&2- 
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After this we finally arrive at 

(2.15) 

Mp{s) 

Mfiis + 1) 

(2.16) 

r(l + + /3‘^{s + a 2 ))r(l + + /3^(<s + 62 ))r(l + /3‘^{s — l))r (2 + fli + + (s + 1 + 02 + ^ 2 )/^^) 

r(l + /3^(s + 1/2 + 02/2 + &2/2) + (oi + 6i)/2)r(/3^(s + (1 + 02 + &2)/2) + (oi + &i + 3)/2) 

^2~2(l+ai +foi) —(4s+l+2a2 +2^2 

r(l + (oi + &i)/2 + f 3 ^{s + (02 + 62)/2))r((3 + Oi + bi)/2 + / 3 ^(s + (02 + b2)/2)) 


To determine the function Mp{s) which satishes fl2.16p for any complex s we follow [25] and 
introduce a variant of the Barnes function Gp{x) which for any 3fJ()r) > 0 is dehned by: 


, ^ , x-Q/2, , , rdt 

lnG'; 3 (x) =- - -ln(27r) + — 


— —t 

n <-) 


-,—xt 


—t 


" _+ £Z(o/2-rf + 

t \ (l-e-0‘)(l-e-W) ^ 2 ’ t 


where Q = (3 + 1/(3. This function satishes the so-called self-duality relation 

(2.17) Gjs^x) = Gijp^x) 

and further posesses a shift property that is central for our studies 

(2.18) Gp{x + (3) = l3^/^-f^^{2Ti)^T{(3x) G^{x) 

One can check that Gj^{x) for (3 = 1 coincides with the standard Barnes function G{x) which 
is a unique solution of the recursion G{x -|- 1) = r(a:)G(a;) satisfying G'(l) = 1 . Similarly to 
the standard Barnes function the general Barnes Gf^^x) has no poles and only zeroes located 
at X = —n(3 — m/(3, 7 T,,m = 0,1,... A detailed discussion of properties of functions closely 
related to Gp{x) can be found in jlOl ITT] . 

Let us now dehne a function (s) of the complex argument s by 


Mf\s) = 


(2.19) 


X 


X 


G /3 {(3 {s + (/3 (s + 

Gfj (^(3{s + 02 ) + Gjs (^(3{s + ^ 2 ) + 

G/3 (/9 (s + (/3 (s + 

G /3 (^( 3 {s -|- 1 -|- 02 -|- ^2) + Qp (^( 3 {s — 1) 


where Bi = 2(3“^ and B 2 = 2(ai -|- 61 -|- 1) - 1 - f3‘^{2a2 -l- 262 — !)• Then a straightforward 
computation which relies on the identity following from (I2.18p 


( 2 . 20 ) 


10 


























shows that the ratio 


( 2 . 21 ) 


Mf\s) 

Mf\s+1) 


reproduces the right-hand side of fl2.16l) from which we conclude 

Mf\s) ^ Mp{s) 

Mf\s + l) Mp{s + l) 


which hnally implies that 

(2.22) Mi,(s) = 

l-Lj 

where = 1. Together with fl2.4p . fl2.5p and the fact that Ma(l) = 

(5<l- 

(2.23) E(Zjv(/3)^~^) - 


1 , we obtain for 

., <T) 

Aff (1) 


2.2. Duality and the freezing transition. The pair (I2.19p - (l2.22p solves the problem of 
finding the complex moments M^(s) = E(z^“®) of the random variable for any complex s, 
and (3 < 1. Knowledge of such moments can be used to restore the probability distribution 
of Zjs, hence of the partition function Z]sf[f3), and of its logarithm (the free energy) for large 
S> 1. Our goal is however to study the limit of the latter as (3 ^ oo and one therefore 
should have a way of extracting information on the distribution for /3 > 1. In doing this 
we rely on the freezing transition seenario for logarithmically correlated random landscapes. 
The background idea of such scenario goes back to [10] and was further advanced and clarified 
in the series of works [211 ESI ESI E2] • In brief, this scenario predicts a phase transition at 
the critical value /9 = 1 and amounts to the following principle: 

Thermodynamic quantities which for (3 <1 are self-dual functions of the inverse 
temperature (3, i.e. functions that remain invariant under the transformation (3 —>• (3~^, 
retain for all (3 > 1 the value they acquired at the point of self-duality f3 = 1. 

Although such a scenario is not yet proven mathematically in full generality and has the 
status of a conjecture supported by physical arguments and available numerics, recently a 
few nontrivial aspects of freezing were verified within rigorous probabilistic analysis, see e.g. 
[SI [151 US] for efforts in this direction. 

Within that scenario, one of the main outcomes of the analysis performed in [25] is that the 
self-dual object associated with the distribution of the partition function for logarithmically 
correlated landscapes is expected to be the appropriately defined Laplace transform: 

(2.24) gf,iy) = E (exp {-e^^yZ^{f3)lZ%{f3)\) , 


where Z^[f3) is a typical scale of the partition function which is extracted from the asymp¬ 
totic for the integer moments and in our case can be chosen as 


(2.25) 




h) 

N 


(B) = W + _(1 

G( 2 / 3 -M)r(l-/ 32 ) 


Moreover, dehning the probability density pg{y) by pyiy) = —g'p{y) one can show that 
the double-sided Laplace transform for such a probability density is related to the complex 
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moments Mp[s) = E 
(see eq.(26) of [25] ) 


Zn{0) 


1 -s 


of the scaled partition function via the following relation 


(2.26) 


/■°° ~ s s 

In J pi 3 {y)e^" dy = \nMis{l +-)+ \nT{l +-) 


Actually, as shown in [25] the freezing scenario implies that the variable y whose probability 
density is given by p^=i(|/) is precisely the fluctuating part of the height of the global min¬ 
imum of the random potential which is our main object of interest. Note however that the 
scale z‘'^\l3) diverges when approaching the critical point /3 = 1 , and that the associated 
free energy —^ \ogZ^^{[d) is self-dual only in the leading order, given by —(/3 -|- /9“^) log A^. 
The latter term after freezing at /3 = 1 yields the leading 21ogA^ term in our conjecture 
Eq. fll.3ll for the maximum, whereas the logarithmically divergent term — ^ logr(l — (3"^) af¬ 
ter careful re-interpretation results in the second term — | log log iV, see |2^ for the detailed 
explanation of that mechanism. The procedure leaves however a certain arbitrariness in the 
terms of the order of unity in the mean free energy, hence in the overall shift of the position 
of the maximum. Let us stress however that apart from such a shift, the shape of the distri¬ 
bution function recovered in the framework of the freezing paradigm is completely hxed by 
the procedure. 

Our strategy therefore will be to check if self-duality holds for the right-hand side combi¬ 
nation in fl2.26p when we substitute our expression for the moments. Before we proceed, it 
will be helpful to further expand our expression fl2.23p . Inserting fl2.19p and making use of 
the identity 


(2.27) 


_\= r(l + (3‘^{S 1)) . ' ,o_i/2-/32(s-1) 


shows that (taking into account all prefactors coming from (12.41) . (12.5p and (I2.19p ) 


E(Z^(/3)1-*) ~ ^ - 


X 



G/3 (/?(^ 

Gy (p (s -f + 3 +y j 


Mf\l) 

's + 22 ±| 2±1 j ^ 2 +aJfei j ^ ^ 

(s -|- 02) + Gy ip{s -\- 62) + 


A direct inspection makes it clear that the self-duality is only possible if either ai = 02 , = &2 

or oi = 1)2, hi = 02 . For the GUE characteristic polynomials, we have Oi = 61 = 1/2, 
02 = = <?/2 so that duality occurs only if q = 1. We therefore have to choose g = 1 to be 

able to rely upon the freezing scenario allowing to interpret the function py=i{y) calculated 
from its Laplace transform via fl2.26p as the probability density for the (shifted) global 
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minimum. Using fl 2 . 28 p with oi = 02 = 6i = 62 = | we get 

(2.29) G*/? (s + 1) + ^ j G;3 (s + 1) + I j G/3 ^/3 (s + I) + I j 

G /3 ^/3(s + 1/2) + ^j ^/l(s + 2) + ^j Gp y(3s + ^j 

where is a constant determined by the condition E(Zjv(/5)^“®)|s=i = 1. Inserting fl2.19p 
into the right-hand side of fl2.26p (which is now manifestly self-dual) leads to the following 
expression at /9 = 1: 

(2.30) 


/ P0=i{y) e^^dy = i^T(l + s)M^=i(l + s) 

J —CO 

(2.31) 

= i_/y*r2n I + 7/2)^G(s -H 3)G(g -F 4) ^ 1 ^^^ ^(l + s)G{s -F 7/2)^r(s -F 3) 

G ^ ’ G2(s + 3)G(s -1- 6)G(s + 2) G G(s + l)G(s + 6) 


where G = c^=i and is a constant which determines the shift in the maximum as discussed 
below fl2.26p . The value K = 271 in fll.ip is conjectured from the results of numerical 
simulations in the next section. The latter formula fl2.3ip constitutes our main analytical 
result and hnally leads to our Conjecture 11.11 


3. Numerical study of the distribution of the maximum modulus of GUE 

CHARACTERISTIC POLYNOMIALS 

The purpose of this Section is to provide a numerical test of Conjecture 11.11 

3.1. Results. In Figure [3T] we present a histogram of the recentered and rescaled maximum 
of the GUE characteristic polynomial, dehned by 

(3.1) y% := (21og(iV) - (3/2) log(log(iV)) - + c^)(l + s%) 

with dehned in fll.2p . Here we used the matrix size N = 3000 and 250, 000 realizations 
of the GUE ensemble. The dashed red line is the exact probability density of the random 
variable y dehned via its Laplace transform in fll.4p . In fl3.ip we have recentered and scaled 
by c*^ = 0.216 and s*^ = 0.188, presumably a consequence of hnite-N ehects due to the o(l) 
terms in fll.3p . Note that the inhuence of shift/recentering is already quite small compared 
with the predicted considerably larger (3/2) log(log(N)) ~ 3.12 shift. The parameters 
and s*^ were calculated empirically from the mean and variance of y in fll.4p according to 
the formula 

4 = ny)l{s% + 1) - (2 log(iV) - (3/2) log(log(iV)) - M%), 

as derived by requiring E(?/^) = E(i/) and Var(i/)(r) = Var(?/). In Table [Uwe display values of 
the parameters c*^ and for the studied range of sizes N, as determined empirically from 
the mean and variance of the random variable u. The observed decay with N is certainly 
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Figure 3.1. The centered and scaled maximum as defined by fl3.1l) . The 
dashed line is the probability density of the random variable y given in Laplace 
space by (na. 
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Figure 3.2. Each triangle represents a value of ^*N obtained from fl3.2p with 
250, 000 realizations. 
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Table 1. Finite-corrections for increasing values of N all with 250,000 realizations 


N 

Cff 

Sn 

150 

0.329 

0.331 

600 

0.267 

0.248 

1050 

0.244 

0.224 

1500 

0.234 

0.212 

1950 

0.228 

0.202 

2400 

0.221 

0.195 

3000 

0.216 

0.188 


0.35 

0.3 

0.25 

0.2 

0.15 

0.1 

0.05 


Figure 3.3. The inverse Laplace transform of formula (11.7p . 

consistent with asymptotic validity of our Conjecture 11.11 though the convergence to the 
asymptotic results is too slow to make more definite claims. To resolve further decrease 
of the coefficients and Sn would require much larger matrices and is computationally 
demanding. 

Finally, we provide a numerical validation of the decomposition fll.hp . In Figure [33] we 
plot the inverse Laplace transform of fll.71) obtained by a direct numerical evaluation of the 
integral in the Bromwich inversion formula for the Laplace transform. The positive and 
normalized curve clearly corresponds to a bona fide probability density of some real random 
variable y'. 

3.2. Numerical method. The numerical evaluation of the maximum value fll.2p may be 
considered quite a non-trivial problem in its own right, for at least two reasons. Firstly, 
the characteristic polynomial Pn{x) having zeros as the eigenvalues of "H, displays 0{N) 
oscillations in the spectral interval [—1,1] with hugely varying peaks heights. This produces 
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considerable clusterings of ‘near-maxima’ which may confuse any naive attempt to find the 
true maximum value. Secondly, the slow changing nature of the correction terms in Conjec¬ 
ture [TTTl of order log(iV) and loglog(iV)) respectively, require one to go to somewhat large 
matrices to resolve reasonable asymptotic behaviour. The problem is further compounded 
by the numerical instability of calculating determinants of such matrices. 

Our solution to these problems heavily relies on a sparse realization of CUE matrices H 
originally due to Trotter mi (see also Dumitriu and Edelman HZ])- He discovered that 
the eigenvalues of CUE matrices H have the same joint probability density as those of the 

\ 

X 4 

X2{N-2) A/’(0,2) X2(N-1) 

X2{N-1) A/’(0,2) / 

where A/'(0, 2) is a normal random variable with mean 0 and variance 2. The sub-diagonal 
is composed of random variables X 2 n having the same density as where xln is a 

y-square random variable with 2n degrees of freedom. To compute the maximum value of 
Pn{x) = det{xl — H) = det{xl — Ti), we begin by exploiting the known asymptotic behaviour 

(3.4) 2E log \pn{x) I = N{2x^ — 1 — 2 log(2)) -h o(l) 
so that 

(3.5) /Ar(x) := 21og|pAr(x)| - 2E log |p 7 v(x) | ~ 21og|det(e"("’""^/^"'°*^(^ii(a:/-'H))| 

Further progress is now possible thanks to the fact that determinants of tri-diagonal matrices 
satisfy a linear recurrence relation. Furthermore, by an appropriate rescaling, the recursion 
computes determinants of all leading principal minors simultaneously, thus computing fj{x) 
for all j = 1,N in linear time. 

Now to find the maximum, we define a mesh A4 = {—1 + n/A : n = 0,..., 2A} with 
A ~ 2N and evaluate fN^x) at each of the points in M.. At those points where /n^x) is 
maximal the Matlab function ‘fminbnd’ is invoked to converge onto the global maximum. 
Figure [T^ illustrates the complexity of the problem. Our algorithm is sufficiently precise to 
distinguish the true maximum (located at x ~ —0.3 in red) from other possible candidates, 
e.g. X ~ —0.7 as well as the thousands of other local maxima. 
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